Tunneling Spectroscopy of Quantum Charge Fluctuations in the Coulomb Blockade 
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We present a theory of Coulomb blockade oscillations in tunneling through a pair of quantum dots 
connected by a tunable tunneling junction. The positions and amplitudes of peaks in the linear 
conductance are directly related, respectively, to the ground state energy and to the dynamics 
of charge fluctuations. We study analytically both strong and weak interdot tunneling. As the 
tunneling decreases, the period of the peaks doubles, as observed experimentally. In the strong 
tunneling limit, we predict a striking power law temperature dependence of the peak amplitudes. 



The charge of an isolated conductor is quantized in 
units of the elementary charge e. Surprisingly, even if 
the conductor is connected to a particle reservoir by 
a tunnel junction, its charge can still be almost quan- 
tized at low temperatures, a phenomenon known as the 
Coulomb blockade [0J|]. The simplest system which 
shows a Coulomb blockade consists of a small metal- 
lic grain separated from a bulk lead by a thin dielec- 
tric layer. An electron tunneling through the layer in- 
evitably charges the grain, thus increasing its energy by 
Ec = e 2 /2C, where C is the capacitance of the grain. At 
temperatures T <C Ec a negligible fraction of the elec- 
trons in the lead have an energy of order Ec, and one 
might expect that no tunneling into the grain is possible. 
More careful consideration shows, however, that even at 
T = the electrons can tunnel to the virtual states in 
the grain, thus lowering the ground state energy of the 
system ||. Due to this virtual tunneling, the average 
charge of the grain is no longer quantized and acquires 
a correction proportional to the conductance of the bar- 
rier. Charge quantization is completely destroyed when 
the conductance of the barrier approaches e 2 /h {|,|U|. 
Unfortunately a direct measurement of the equilibrium 
properties, such as the average grain charge, comprises a 
challenging, though not impossibl e |?|, e xperiment. 

Several recent experiments |||],[L0] ll]|l^] have probed 
some equilibrium properties of a Coulomb blockade sys- 
tem by measuring the tunneling conductance through a 
pair of coupled quantum dots. Focusing on the exper- 
iment of Waugh, et al. ||, we develop in this paper a 
quantitative theory of the linear conductance in such a 
system; in particular, we predict the gate voltage and 
temperature dependence of the conductance. 

We will study the properties of a two-dot system which 
is schematically shown in Fig. The electrostatic energy 
of this system is a quadratic form of three variables: the 
charges of each dot, eN\ and eN 2 , and the gate voltage 
V g . In the most general case, this energy can be written 
in the following form: 



Here X is a dimensionless parameter proportional to V g , 
and the constants Ec, Ec, A, and a are determined by 
the geometry of the system Q . In Eq. ([!]) we expressed 
the energy in terms of the total number of particles in 
the two dots Ni + N 2 and the relative number Ni — 
N 2 ■ These variables are convenient because the former is 
constant in the absence of tunneling into the leads, and 
the latter describes charge fluctuations between the dots. 
In this paper we will concentrate on the case of symmetric 
geometry of the system, corresponding to A = a — 0, 
which is apparently the case in the experiment Q . 

We first discuss the location of the peaks in the con- 
ductance G of the double dot system. In the limit of very 
small inter-dot conductance, Go <C e 2 /h, the peaks in G 
occur when the electrostatic energy is degenerate; that is, 
when U(n + 1, n) equals either U(n, n) or U(n + 1, n + 1) 
where n is an arbitrary integer. As a result we find peaks 
at the following sequence of gate voltages: 




(2) 



Weak electrostatic coupling between the dots corre- 
sponds to Ec — Ec <C Ec- In this case the two peaks 
with the same n in sequence (^|) merge. This limit is 
observed in the experiment 

If all the junctions shown in Fig. [I] have small conduc- 
tances, then the charges of both grains are well defined. 
In order to study the effect of quantum fluctuations on 
the ground state energy of the system, it is enough to 
increase only the inter-dot conductance Go, keeping all 
other conductances small. Under these conditions, the 
sum J\f = N\ + N2 is constant and can still be treated as 
a c- number, but Ni — N 2 starts to fluctuate. These fluctu- 
ations change the ground state energy, denoted £V, from 
the electrostatic estimate. The peaks in G now occur at 
gate voltages X where Ejy(X) = Ej^ + i{X). In order to 
find Ej^f, one should consider the quantum mechanical 
problem with the Hamiltonian 
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Here the terms Hq and Ht describe, respectively, free 
electrons in the dots and tunneling between the dots; we 
have replaced N2 by M — N\ , and Ni should be treated 
as a quantum operator N± . Typically the size of the dots 
exceeds the effective Bohr radius (~ 100A for GaAs), 
and therefore the level spacing for electrons in the dots is 
much smaller than the charging energy. We will neglect 
the level spacing and assume a continuous spectrum in 
H , in contrast to Refs. [jl4|,[l5|. In the continuous model, 
the non-interacting part of the Hamiltonian, Hq + Ht, 
does not depend on the total number of particles A/". 

In the Hamiltonian (^) the parameter M is an integer. 
However, formally we can consider Eq. @ at any N ' , 
enabling us to relate the ground state energy Ej*/ to the 
average value Ni(ftf) of the first dot's charge: 



dE M 
dN 



2E C (M - 2X) - AE C 
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The condition for the peak position, -Ejv+i — E^ = 0, can 
now be obtained by integration of Eq. (El) . The result is 
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where now Af is an integer again. In the limit Go — > 0, the 
average Ni(Af') is the integer closest to Af'/2, and taking 
M =2n and M = 2n + 1 in Eq. (||) reproduces Eq. @). 
The advantage of Eq. (||) is that it is valid at any Go- 
However, to make use of it, one has to evaluate Ni(Af) 
at any A/", which is a challenging quantum mechanical 
problem because of the Coulomb interaction in Eq. (|^). 
Fortunately, in the limit of a continuous spectrum, the 
Hamiltonian (|3|) coincides with the one for a single dot 
connected by a tunnel junction to a massive lead. The 
latter problem has been extensively studied in the limits 
of weak and strong tunneling into the dot P,^,p|,|6|Jl6[ . 

For weak tunneling, Go <C e 2 /fi, the deviation of 
N±(J\f) from an integer is small and can be found || 
from second order perturbation theory in Ht- For J\f in 
the interval (2n — 1, 2n + 1), we get 



hG , M-2n + l 
iVi(A0 = n + — -j In - r . 



(6) 



Substituting Eq. (^) into Eq. (||) and using the periodicity 
of Ni(Af), we find peak positions shifted by tunneling: 
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The splitting of the two peaks with the same n grows 
linearly with Go- The splitting here results from quan- 
tum charge fluctuations between the dots, not changes in 
geometric capacitances as inferred in Refs. ]lC|Jll|] . 

Clearly, the charge fluctuations grow with Go. In the 
limit of strong tunneling the discreteness of charge N\ 



is no longer important H||@], and Ni(J\f) — » \M [see 
Eq. (0)]. As a result the peaks are equidistant, X* = 
(2A"+ l)/4, which is expected because in this limit the 
two dots form a single conductor. The doubling of the 
period of the peaks from this regime to that of Eq. ^ is 
one of the main observations of the experiment of Waugh, 
et al. @. 

To find the peak positions as the system approaches 
the strong-tunneling limit, one has to specify a model 
of the junction between the dots. For electrostatically 
controlled dots in semiconductor heterostructures, the 
junction is a microconstriction with smooth boundaries 
JItJ . The ground state energy of such a system near the 
strong-tunneling limit, when the reflection coefficient for 
the single transverse mode propagating through the con- 
striction is small 72. = 1 — nhGo/e 2 <§; 1, was found in 
Rcf. D [Eq. (48)]. Using that result, we get 
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2M- 
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where C « 0.5772 is Euler's constant. 

Fig. H compares our results (0) and (||) to the observa- 
tions of Waugh, et al. || . Because the different gates are 
not independent, the relation between the gate voltage Vq 
at which the inter-dot conductance is measured and the 
Vo at which the splitting is measured is not known. We 
take this relation to be a rigid shift || and have used the 
data of Ref. g with a shift of -0.005 V. The agreement 
between theory and experiment is very good indeed. 

As we have seen, the positions of peaks in the linear 
conductance carry information only about the ground- 
state energy of the two-dot system. To study the ex- 
citations, one can analyze the heights and shapes of the 
peaks. If the inter-dot tunneling is weak, Go <C e 2 /U, the 
excitation spectrum consists of two independent quasi- 
particle spectra of the two dots. This enables us to apply 
the standard master-equation technique [fi| and find 



G 



GlG r 



4E C (X - X*)/T 



2(Gi + G r ) sinh[4£ c (A - X*)/T] ' 
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Here A* is the position of the center of a peak given by 
one of the values in the sequence (0). In deriving Eq. (|) 
we assumed that the tunneling into the leads is much 
weaker than between the dots, G; ~ G r <C Go- Thus 
Eq. (||) reproduces the result for a single dot [^8|. This 
simple formula is valid only at sufficiently low temper- 
atures, when the width of a peak SX ~ T/Eq is much 
smaller than the spacing between the peaks. For small 
capacitive coupling between the dots {Ec ~ Ec), this 
yields T< (hG a /e 2 )E c . 

In contrast to weak tunneling, at Go ~ e 2 /% the excita- 
tion spectra of the two dots are not independent. In this 
regime an electron tunneling into the left dot shakes up 
the quantum state of the whole two-dot system, leading 
to a suppression of the conductance at low temperature. 
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To illustrate this phenomenon, we calculate the tempera- 
ture dependence of the peak heights in the case of perfect 
inter-dot transmission, Go = e 2 /n%. 

As we have seen, the conductance peaks in the strong- 
tunneling limit are equally spaced, as if the double dot 
system were a single dot. In fact, this remains true even 
for asymmetric double dots [A, a ^ in Eq. ([!])] because 
when 1Z — > the energy is simply given by the first 
term in Eq. ([!]). However, the specific geometry of the 
system — whether it is one or two dots, and the degree of 
asymmetry — will show up in the peak heights. Unlike in 
a single dot, the single-mode constriction impedes charge 
propagation between the two dots, thus producing effects 
similar to those for a single junction coupled to an en- 
vironment When an electron tunnels from the lead 
into the left dot, the other electrons in both dots must re- 
distribute in order to minimize the electrostatic energy: 
a charge of (1 + A)/2 electrons must pass through the 
constriction. As a result, the overlap of the two ground 
states, before and after the tunneling, vanishes, as in An- 
derson's orthogonality catastrophe. 

At non-zero temperature, the tunneling density of 
states is suppressed as T 7 , where the exponent is re- 
lated to the scattering phase shifts 5 m in each one- 
dimensional channel m by 7 = ^2 m (8 m /ir) 2 . Accord- 
ing to the Friedel sum rule, 6 m /ir is the average charge 
transferred into each channel. A single-mode constriction 
provides two channels for each dot (two spins), yielding 
4 channels in total. In our case, S m /Tr = ±(1 + A)/4, 
where the plus (minus) sign is for the channels in the 
right (left) dot. Thus, the rate of tunneling into the left 
dot is suppressed by the factor T^ 1+x ^ / 4 . For the rate 
between the right dot and lead, one should replace A by 
—A. Because the junctions are connected in series, the 
smaller of the two rates determines the conductance, 



G = G h 
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(10) 



Here the coefficient Gb is of order G; ~ G r , and the 
peaks are centered at X* = (2AT+ 1)/4. For a symmetric 
system, A — 0, the peak conductance obeys G oc T 1 / 4 . 
The temperature dependence in Eq. (nm can be obtained 



analytically 20 in the spirit of the bosonization approach 
p9|,p[; this technique also yields for the peak shapes 



F 7 (;r) 



1 l r (i + i + fr) 

cosh(x/2) T(2 + 7) 
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At 7 = the shape given by Eq. (|llj) is identical to 
that for weak tunneling Eq. (|J): in both cases there is 
no charge transfer between the dots during the act of 
tunneling, and the spectra of the two dots are decoupled. 

Comparing our results in the weak and strong tunnel- 
ing limits [Eqs. (||) and (M)], we see that when Go grows 



the conductance of the system decreases due to the or- 
thogonality catastrophe. For inter-dot conductances be- 
tween these limits, the system must crossover from tem- 
perature independent peak heights for weak tunneling to 
the power-law suppression of the conductance for strong 
tunneling. The theory of this crossover will be reported 
elsewhere pf| . 

We have seen that at Go = e 2 /irk an asymmetry of 
capacitances, A ^ 0, causes temperature dependent peak 
heights (|l0|). In this limit, the peaks are equidistant and 
have the same height. We will now show that for weak 
tunneling even a small asymmetry dramatically affects 
the whole pattern of peaks. Below we study the asym- 
metry related to a small non-zero a in the electrostatic 
energy ([!]) ; the results are easily generalized to include a 
small A by replacing a — > a + 2A. Motivated by experi- 
ment [||], we will assume Ec — Ec- 

First, let us determine the positions of the peaks in 
conductance. The condition for a peak is the degeneracy 
of the energy (|l|) with respect to adding an electron to 
either of the dots; this yields the two sequences of peaks 



X* = 



1/2 



l-a/2' 



XX = 



1/2 



l + a/2 



(12) 



Since a <C 1, the two periods are very close, and one 
observes beats with an approximate superperiod of a . 

The asymmetry lifts the degeneracy of the two states 
with an extra electron on either the left or right 
dot. The energy gap between these two states is 
A(A) = 4Ec\aX — m\, where m is the integer nearest 
to aX. If the temperature is much lower that A, tun- 
neling between two real states is suppressed. Instead, an 
extra electron in the left dot can escape into the right 
lead through a virtual state in the right dot. Because we 
assume that the level spacing is small compared to tem- 
perature, the dominant escape mechanism is inelastic co- 
tunneling pl[|, i.e., an extra electron-hole pair is created 
in the right dot. The total rate |2l| of such processes 
is proportional to (T/A) 2 and limits the conductance at 
sufficiently low temperatures. The calculation of peaks 
in conductance |20| leads to 



G 



T 



A(X) 



AE C (X - X*] 
T 



(13) 



Here the peak positions X* are defined by Eq. (Jl 2j) . The 
oscillations of the energy gap A(X) modulate the peak 
heights with the period oT x in gate voltage X. 

In conclusion, the fluctuation of electron charge be- 
tween two quantum dots strongly affects tunneling 
through such a structure. First, the conductance peaks 
are split because of the lowering of the ground state en- 
ergy by charge fluctuations. This splitting, and eventual 
halving of the period of the conductance peaks as the 
interdot conductance grows, is a dramatic feature of the 
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experimental data H which is fully confirmed in our the- 
ory. Second, the temperature dependence of the peak 
height and shape is directly related to the dynamics of 
the quantum charge fluctuations. For a double-dot con- 
nected by a reflectionless constriction, this produces a 
striking fractional power law temperature dependence. 
Our theory is valid for a wide range of temperatures, 
limited only by the charging energy from above, and by 
the discrete energy level spacing from below. 
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FIG. 2. The normalized splitting of the Coulomb blockade 
peaks as a function of the inter-dot conductance. Our theo- 
retical results (dashed lines) are in good agreement with the 
experiment of Ref. ^| [x for data of Fig. 3(a) and + for Fig. 
3(b)]. The splitting is normalized by the period of the peaks 
in the strong tunneling limit. A small asymmetry, a = 0.05, 
has been included based on the experimental parameters. 
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Figure 2. 



